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COMPUTATIONAL  METHODS  FOR  BIFURCATION  PROBLEMS 
WITH  SYMMETRIES  ON  THE  MANIFOLD* 

BIN  HONGt 

Abstract  :  This  paper  is  about  numerical  methods  for  the  determina¬ 
tion  of  bifurcation  points  of  certain  steady  state  multi-parameter  prob¬ 
lems  in  the  presence  of  symmetries.  A  principal  tool  is  the  fact  that 
under  general  conditions  the  solution  set  forms  a  manifold  in  the  space 
of  all  state  and  parameter  variables.  The  reduced  manifold  with  respect 
to  some  subsymmetry  is  introduced.  Methods  are  presented  for  the  lo¬ 
cal  computation  of  the  submanifold  of  bifurcation  points  of  the  same 
symmetry. 

Key  words  :  solution  manifold,  symmetry  groups,  reduced  form,  fold- 
point  calculation 

AMS(MOS)  subject  classification  :  65J15,  58F14 

1.  Introduction.  Equilibrium  problems  for  many  physical  systems  are  modeled  by 
multi-parameter  dependent  nonlinear  equations 

F(z,  A)  =  0.  (1.1) 

Here  z  varies  in  some  space  Z  and  characterizes  the  state  of  the  system  while  A  denotes  the 
parameter  variable  allowed  to  vary  in  a  space  A  of  dimension  p  >  1.  Under  fairly  general 
conditions,  the  regular  solution  set  of  (1.1)  forms  a  p-dimensional  differentiable  manifold 
and  interest  centers  on  determining  the  foldpoints  of  the  manifold  M  with  respect  to  the 
parameter  space  A  (see  [19]). 

Often  in  physical  applications  the  system  (1.1)  is  covariant  with  respect  to  a  trans¬ 
formation  group  G;  that  is, 

F(Tz(g)z ,  Tx(g)\))  =  S(g)F(z,  A),  for  all  g  €  C,  (1.2) 

where  T-,  T\  and  S  are  group  representations  of  G.  When  such  a  symmetry  is  present, 
it  can  aid  considerably  in  the  computation  of  determination  of  the  bifurcation  structure, 
especially  at  multiple  bifurcation  points.  The  computation  of  a  multiple  bifurcation  point, 
in  general,  is  difficult  and  costly.  A  major  advantage  of  applying  group  theoretic  methods 
is  the  possibility  of  considering  the  problem  in  a  reduced  form  reflecting  some  subsymmetry 
under  which  certain  multiple  bifurcation  points  reduce  to  the  simple  case  and  hence  can 
be  computed  by  algorithms  available  for  such  simple  bifurcation  points. 

*  This  work  was  in  part  supported  by  the  National  Science  Foundation  under  grant 
CCR-8907654  and  the  Office  of  Naval  Research  under  grant  ONR-N-00014-90-JT025. 

f  Institute  for  Computational  Mathematics  and  Applications,  University  of  Pittsburgh, 
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In  recent  years,  a  number  of  papers  considered  symmetry  and  symmetry  breaking  in 
connection  with  the  computation  of  bifurcations  (see  e.g.  [5],  [12]).  But  most  of  them 
are  restricted  to  one-parameter  problems,  and  even  those  concerned  with  multi-parameter 
problems  single  out  a  specific  parameter  for  the  computation  and  hold  all  other  parameters 
at  a  fixed  value  (see  e.g.  [3]).  The  aim  of  this  presentation  is  to  use  group-theoretic  tech¬ 
niques  together  with  a  differential-geometric  viewpoint  to  analyse  bifurcation  phenomena 
of  multi-parameter  problems  without  giving  any  preference  to  any  one  parameter. 

Section  2  summarizes  the  needed  background  material.  We  adopt  from  [19]  and  [7] 
the  differential-geometric  setting  for  the  solution  manifolds  and  then  introduce  some  basic 
concepts  of  group  theory.  This  includes  some  material  on  the  discretization  of  infinite¬ 
dimensional  problems  under  symmetry. 

Our  approach  is  based  on  previous  work  in  the  setting  of  solution  manifolds  such  as 
methods  for  following  a  path  on  M  (see  [19])  and  for  computing  simplicial  approximations 
of  subsets  of  M  (see  [20]),  the  geometric  framework  for  the  numerical  study  of  singular 
points  introduced  in  [7]  and  a  resulting  method  for  the  computation  of  some  type  of  singular 
points  (see  [4]).  The  key  objective  of  this  research  is  to  incorporate  group  actions  into  these 
numerical  methods. 

Section  3  concerns  group  actions  on  solution  manifolds.  In  Section  3,  the  reduced 
manifold  M //  with  respect  to  some  subsymmetry  H  is  introduced  and  the  relation  between 
the  original  solution  manifold  M  and  the  reduced  solution  manifold  Mh  is  discussed.  It  is 
shown  that  the  usual  trivial  group  action  on  the  parameter  space  is  not  always  applicable 
here  and  that  a  more  general  one  has  to  be  introduced.  For  this,  group  actions  on  the 
parameter  space  A  are  considered  in  the  case  of  standard  unfoldings  in  bifurcation  theory. 
For  the  study  of  symmetry  breaking  bifurcations  the  group  action  is  extended  to  the  so- 
called  cut  function  and  its  (k  - f  l)-form  which  in  [7]  form  the  basis  of  the  framework  for 
the  study  of  bifurcation  phenomena  on  manifolds.  An  important  question  in  the  study  of 
symmetry  breaking  is  the  determination  of  the  choice  of  the  symmetry  breaking  subgroup. 
For  this,  we  present  a  manifold  version  of  the  equivariant  branching  lemma  which  provides 
some  choices  for  the  symmetry  breaking  subgroup.  In  other  work  on  the  application  of 
group  representation  theory  to  the  computation  of  bifurcations  (see  e.g.  [5],  [12]),  the 
symmetry  breaking  subgroup  has  to  be  pre-selected  for  which  the  bifurcation  point  is 
to  be  found.  We  discuss  here  how  to  implement  numerically  the  equivariant  branching 
lemma  to  provide  better  information  for  detecting  the  symmetry  breaking  subgroup  at 
a  bifurcation  point.  A  recent  theory  about  the  selection  of  the  choice  of  the  symmetry 
breaking  subgroups  developed  by  Rabier  ([17],  [18])  is  shown  to  be  of  considerable  help 
as  well.  By  means  of  the  equivariant  branching  lemma;  that  is,  by  consideration  of  the 
original  problem  in  a  reduced  form  reflecting  some  subsymmetry  H ,  it  becomes  possible 
to  reduce  certain  multiple  singular  points  to  simple  ones  and  to  apply  the  method  for  the 
computation  of  such  simple  singular  points  given  in  [4].  The  set  of  these  singular  points 
with  the  same  symmetry  H  turns  out  to  form  a  submanifold  of  the  reduced  manifold  M//. 

Finally  in  Section  4,  we  illustrate  the  behavior  of  the  methods  of  Section  3  with  an 
example.  In  this  example  comparison  of  the  reduced  manifold  and  of  the  submanifolds  of 
bifurcation  points  at  a  multiple  symmetry  breaking  bifurcation  point  reflects  a  rich  bifur¬ 
cation  behavior  and  indicates  also  some  interesting  facts  about  nonconjugate  symmetry 
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breaking  bifurcation  points. 


2.  Background.  Throughout  the  paper  the  following  assumption  shall  hold. 

Assumption  2.1  s  (i)  X,  Y  are  real  Hilbert  spaces;  (ii)  D  is  an  open  subset  of  X;  (iii) 
F  \  D  Y  is  a  Fredholm  mapping  of  class  Cr,r>  1  and  index  p  >  1. 


We  consider  the  equation 


F(x)  =  0, 


(2.1) 


and  its  regular  solution  set 


M  =  {i:xE  S,F(x)  =  0,rge  DF(x )  =  F}, 


(2.2) 


which  will  always  be  assumed  to  be  nonempty.  It  is  well  known  that  M  is  a  p-dimensional 
Cr-manifold  in  X  without  boundary.  For  simplicity,  the  tangent  space  TZoM  at  the  point 
xo  €  M  of  this  manifold  will  be  identified  with  the  kernel  of  the  derivative  DF(x o)  of  F; 
that  is,  we  set 

TXoM  =  ker  DF(x0)  =  {u  :  u  e  X,  DF(x0)u  =  0},  (2.3) 

whence 

NX0M  =  {TX,M)X  =  (ker -DF(xo))'L  =  rge  L>F(x0)*  (2.4) 

may  be  called  the  normal  space  at  the  point. 

A  given  p-dimensional  subspace  W  of  X  induces  a  local  coordinate  system  of  M  at 

xo  6  M  if 

w  n  NXqM  =  {0}.  (2.5) 

When  (2.5)  holds,  it  is  always  possible  to  define  an  isomorphism  A  :  Y  IF_L  of  F  onto 
IF-*-.  Then  there  exists  an  open  ball  J  =  B{ 0,  (3)  C  W  with  radius  0  >  0,  an  open 
neighborhood  U  C  X  of  xo,  and  a  unique  Cr-function  rj  :  J  i— ►  F  such  that  77(0)  =  0  and 

M  HU  =  {x€Ar;x  =  x(t)  =  xq  +  t  +  Ap(t),  Vt  6  J},  (2.6) 


(see  e.g.  [19]). 

As  indicated  already  by  the  form  of  (1.1),  many  applications  involve  a  natural  orthog¬ 
onal  splitting 

A’  =  Z  ®  A,  Z  =  Ax,  dim  A  =  p  (2-7) 

of  the  domain  space  X  into  a  state  subspace  Z  and  a  parameter  subspace  A.  Then  interest 
often  centers  on  determining  the  singular  points  with  respect  to  A.  Let 


Zq  —  Z  n  TXo^l,  No  —  An  NXoM,  (2-8) 

then,  a  point  x0  €  M  is  a  foldpoint  with  respect  to  A’  =  Z  ©  A  if  the  index  q  =  dim  Z0  = 
dim  No  is  nonzero.  This  integer  q  is  the  first  singularity  index  of  x<)  (see  [7]). 
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Let  GL(X)  be  the  group  of  all  nonsingular  linear  transformations  of  X  onto  itself.  A 
representation  of  a  group  G  on  a  vector  space  X  is  a  homomorphism  T  :  g  i— ►  T(g)  of  G 
into  GL(X),  such  that 

T(9l92)  =  T(gi)T(g2),Vgug2  €  G.  (2.9) 

In  this  paper,  we  consider  only  finite  or  compact  groups  G.  Note  that  for  the  group 
representation  T  of  finite  or  compact  groups  G  on  a  Hilbert  space  AT,  it  is  always  possible 
to  assume  T  to  be  unitary,  since  T  is  equivalent  to  a  unitary  representation  (see  [21]). 

The  natural  orthogonal  splitting  (2.7)  provides  us  with  an  isomorphism  A  :  Y  >— ►  Z 
of  Y  onto  Z.  This  isomorphism  of  Y  onto  Z  again  denoted  by  A  in  turn  induces  a 
representation  5  of  G  on  Y  namely 

S(g)  =  A~'T(g)A,Vg  €  G.  (2.10) 

This  requires  that 

T(g)ZcZ,VgeG,  (2.11) 

which  is  very  natural  in  applications.  We  shall  assume  that  S  is  unitary  as  well. 

In  line  with  this  we  say  that  the  mapping  F  of  Assumption  2.1  is  equivariant  under  the 
group  G  or  that  a  covariant  group  action  is  defined  for  the  equation  (2.1)  if  the  following 
four  conditions  hold: 

Assumption  2.2  :  (i)  G  is  a  finite  or  compact  group;  (ii)  T  is  a  unitary  group  represen¬ 
tation  of  G  on  X  such  that  (2.11)  holds;  (iii)  the  representation  S  of  (2.10)  induced  by  T 
on  Y  is  unitary;  (iv)  the  relation 

F(T(j)x)=%)F(i),V?6G.  (2.12) 


holds. 

As  usual,  the  orbit  of  a  point  aro  €  X  under  a  group  G  is  the  set  of  points 

Oxo  =  {x  €  A  :  x  =  T(g)x0yg  €  G}.  (2.13) 

For  any  point  x  6  X  it  is  wek  known  that  the  set 

Gz  =  {g  €  G  :  T(g)x  =  x}  (2.14) 

is  a  subgroup  of  G.  Any  subgroup  H  of  G  such  that  H  =  Gx  for  some  a:  6  AT  is  called  an 
isotropy  subgroup  of  G.  The  isotropy  subgroups  of  the  points  of  any  orbit  are  conjugate 
to  each  other,  that  is 

Gr(g)z  =  gGzg  1,g€G,x£  X.  (215) 

For  any  subgroup  H  C  G,  the  set 

XH  =  {x  €  X  :  T(h)x  =  x,\fh  €  H}.  (2.16) 

is  a  subspace  of  X  called  the  fixed  point  subspace  of  H  in  AT.  We  can  analogously  define 
the  fixed-point  subspace  of  H  in  Y.  Note  that  Hi  D  H2  implies  that  XH'  C  XH 2.  From 
(2.12)  it  follows  that  F( XH)  C  YH .  This  leads  to  the  following  important  conclusion  : 
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Theorem  2.3  :  Let  Fh  denote  the  restriction  of  F  to  XH .  Then  Xq  €  X H  is  a  solution 
of  (2.1)  if  and  only  if  it  is  a  solution  of  the  H -reduced  problem 

Fh(x)  =  0.  (2.17) 

If  X  is  a  finite-dimensional  space  then  (2.17)  is  a  lower-dimensional  system  of  equa¬ 
tions  than  (2.1).  For  the  computation  this  is  a  most  compelling  consequence  of  Theorem 
2.3.  Some  information  about  the  dimension  of  XH  and  the  orthogonal  projection  from  X 
onto  XH  is  given  by  the  following  result. 

Theorem  2.4  :  Let  G  be  a  finite  or  compact  group  acting  on  X  and  let  H  C  G  be  a 
finite  subgroup.  Then  we  have 

£  tr  (T(h)).  (2.18) 

1  1  hen 

Furthermore,  the  linear  operator  Ph  on  X  defined  by 

p«  =  IVT  E  T(h)’  (2'19) 

1  1  heH 

is  the  orthogonal  projection  from  X  onto  XH . 

Theorem  2.4  can  be  generalized  to  compact  subgroups  H.  Then  the  sums  in  (2.18), 
(2.19)  over  the  elements  of  H  have  to  be  replaced  by  Haar  integrals  over  H.  For  this 
generalization  and  the  proof  of  Theorem  2.4  we  refer  to  [11]. 

In  practice,  the  parameter  space  is  always  finite  dimensional,  but  the  state-space  may 
be  infinite-dimensional.  In  that  case,  for  the  computation,  the  mapping  F  in  (2.1)  needs 
to  be  approximated  by  a  finite-dimensional  analogue.  For  this,  we  follow  the  presentation 
in  [6]  and  choose  a  finite-dimensional  subspace  Z/,  of  Z,  a  discretization  projection  Ph  6 
L(Z)  such  that  rg ePh  =  Z/,.  Here  h  belongs  to  some  index  set  of  small  positive  scalars. 
With  Xh  =  Zh  ®  A,Ffc  =  A~*Zh  C  Y ,  and  the  projection  Qh  :  Y  ►  Yh  specified  by 
Qh  =  A~lPhA ,  the  discretized  function  Fh  from  the  open  subset  D  into  Yh  is  then  defined 

by 

Fh  :  D  *  Yh,Fh(x)  =  QhF(x),x  e  D.  (2.20) 

For  Fh  to  preserve  the  equi variance  of  F  under  the  group  G,  we  need  to  know  under  what 
condition  the  group  representation  5  does  commute  with  the  projection  Qh- 

Theorem  2.5  :  ([14,  Theorem  3.11])  If  the  space  X  can  be  decomposed  in  the  form 
X  =  W\  ®  W2,  where  W\  and  W2  are  invariant  under  the  representation  T  of  the  group 
G,  then  the  projection  operator  P  on  X  defined  by 


Pv  =  uij,  for  v  =  uq  +  tU2,Wi  €  , 


(2.21) 


satisfies 

T{g)P  =  PT(g),Vg  €  G.  (2.22) 


5 


Conversely,  if  P  is  a  projection  operator  on  X  satisfying  (2.21 )  then  X  —  Wx  0  W2,  where 
W\  =  rge  P  and  W2  =  ker  P  are  invariant  under  T. 

By  Theorem  2.5,  if  Zh  is  invariant  under  the  group  representation  T  of  G  on  X,  then  T 
is  interchangeable  with  P*.  Moreover,  inherits  from  Zh  the  property  of  being  invariant 
under  the  induced  representation  5  of  G  on  Y.  Thus  S(g)  commutes  with  Qh  as  well. 
With  the  finite  dimensional  representations  Th{g )  =  T(y)|xh ,  Sh(g)  =  *S'(^)| v*,,  of  G  on  Xh 
and  Yh,  respectively,  we  have  the  ’’discretized”  group  action 

Fh(Th(g)x )  =  Sh(g)Fh( x),Vx  €  Xh,Vg  G  G,  (2.23) 

3.  Reduced  manifolds  and  symmetry  bifurcation.  We  assume  in  this  section  that 
our  problem  has  been  brought  into  finite  dimensional  form 

F(x)  =  0,  F  :  Rn  t-+  Rm,  n  =  m  +  p,p  >  1,  (3-1) 

and  consider  the  derivation  of  the  reduced  manifold  with  respect  to  a  given  isotropy  sub¬ 
group.  By  the  orthogonality  of  the  group  representation  T  we  have 

Tt(g)  =  T(g)-1  =  T(g-1),\fg€G.  (3.2) 

Moreover,  for  all  g  G  G  and  any  x  £  M  it  follows  from  (2.12)  that 

DF(T(g)x)T(g)  =  S(g)DF(x), 

T(gy  DF(T(g)xY  =  DF(xyS(g)\ 

and  by  (3.2)  that 

T(g)DF(T(g-')Iy  =  DF(x)<S(g). 

From  (3.3)  we  obtain  that  M  is  invariant  under  T;  that  is,  T(g)(M )  C  M,\/g  6  G. 

As  pointed  out  in  Section  2  when  F  admits  a  symmetry  behavior  then  it  is  compu¬ 
tationally  advantageous  to  work  with  the  reduced  problem  with  respect  to  some  isotropy 
subgroup  H  of  G.  Let  dim(XH)  =  n//  <  dim  A",  dim(FH)  =  m//,  and  Fh  :  XH  Y ** 
the  restriction  (2.17)  of  F  to  XG .  Then  we  define  the  corresponding  reduced  manifold  as 
the  regular  solution  set  of  F//;  that  is, 

Mu  =  {x  :  x  €  XH ,FH(x)  =  0,rge  DFh(x)  -  Y11}.  (3.5) 

Clearly  its  dimension  is  pu  =  nu  —  rnu  and  the  following  theorem  provides  a  connection 
between  Mu  and  the  original  manifold  M. 

Theorem  3.1  :  If  x  6  M  n  XH  then  x  6  Mu- 

Proof:  We  need  to  show  that  for  any  y  €  YH  there  exists  a  v  6  XH  such  that  DFu{x)v  = 
y.  Because  of  x  €  M  we  know  that  for  any  y  G  YH  there  exists  a  Dj  6  X  for  which 
DFu(x)v  1  =  y.  Since,  in  addition,  x  G  XH ,  it  also  follows  that 

DF(x)T{h)v j  =  S(h)y,'ih  G  H.  (3.6) 

By  applying  (3.6)  to  every  group  element  h  G  H ,  and  summing  the  terms  we  obtain 
DF(x)Puv  1  =  Quy  =  V  and  v  =  Puvi  €  XH  which  proves  the  theorem.  | 

The  following  result  is  now  evident  : 


(3.3) 

(3.4) 
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Lemma  3.2  :  For  each  x  €  XH ,  the  null  spaces  ker  DF(x)  and  ker  DF(x)1  are  invariant 
under  T(h )  and  S(h),  respectively,  for  each  h  €  H . 

For  the  study  of  the  solution  manifold  M  it  became  necessary  to  introduce  a  group 
action  on  the  parameter  space  A.  The  usual  trivial  action  on  the  parameter  space,  used, 
e.g.,  in  [10],  [11],  is  not  always  applicable  here  as  the  example  of  the  cusp  function 

F(z,  A,  p)  =  z3  -  Xz  +  p,  (z,  A,  p)T  €  R3,  (3.7) 

shows  which  is  covariant  under  the  ^-symmetry, 

F(-z,  A,  -p)  =  - F(z ,  A,  p).  (3.8) 

It  is  useful  to  consider  the  group  action  on  the  parameter  space  A  in  the  case  of  unfoldings 
of  the  type  standard  in  bifurcation  theory. 

Suppose  that  F  :  Rm  x  R1  Rm  is  a  Cr-mapping  on  some  open  subset  D  containing 
a  point  (20,^0)  €  Rm  x  R1  where  F(z0,  po)  =  0  and  dim  ker  DF(z0,  po)  =  r-f  1  >  2.  These 
conditions  imply  that  dimrge  DF(z0,p0 )  =  m—r  <  m—  1.  Let  a1? . . .  ,ar  6  Rm  be  linearly 
independent  vectors  which  span  ker  DF(z0,  po)1  and  hence  which  span  a  complement  of 
rge  DF{zq,  po)  in  Rm.  Then  an  unfolding  of  F  is  given  by 

F  :  R"  =  Rm  x  R1  x  Rr  ^  Km,F{z,p,6)  =  F(z,p)  +  6iai  +  . . .  +  6rar  (3.9) 

where  6  =  (61, . . .  ,8r)  and  A  =  R1  x  Rr.  By  construction  we  have  rank£>F(a:o,  A0)  =  m 
and  evidently  with  n  =  m  +  1  +  r,p  =  r  +  1,  A  =  (p,  6),  x  =  (z,  A),  the  Assumption  2.1  is 
satisfied. 

Suppose  further  that  the  original  mapping  F  is  equivariant  with  respect  to  a  group 
G ;  that  is,  that 


F(Tz{g)z,p)  =  S(g)F(z,p),Vg  eG,ze  Rm,^  6  R1,  (3.10) 

where  TZ,S  are  group  representations  on  the  spaces  Z  and  Y  respectively.  Moreover, 
assume  that  (zo,p0)  belongs  to  the  fixed-point  space  of  G  on  Rm  x  R1.  We  wish  to  ensure 
that  the  unfolded  mapping  F  of  (3.9)  has  the  same  property;  that  is,  that 

F(Tz(g)z,Tx(g)X)  =  S(g)F(z,  A),V<?  €  G,z  €  Rm,A  €  Rp,  (3.11) 

where  T\  is  a  group  representation  on  the  space  A.  In  order  to  determine  T\(g)  for  any 
g  €  G,  we  multiply  both  sides  of  (3.9)  by  S(g), 

S(g)F(z ,  A)  =  S(g)F(z,p)  +  +  . . .  +  SrS(g)ar.  (3.12) 

By  Lemma  3.2  we  have  S(g)a,  6  span(aj, . . . ,  ar),  so  that 

S(g)at  =  bn(g)ai  +  bi2{g)a2  +  ...  +  bir(g)ar,i  =  l,...,r.  (3.13) 
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Then  the  matrix  B(g )  =  ( bij(g ))  satisfies  det  B(g)  ^  0.  Hence  we  can  rewrite  (3.12)  as 


S(g)F(z,  A)  =  S(g)F(z,  p)  +  S[ax  +  . . .  +  6*rar 

where 


fs'A  /s‘\ 


By  combining  (3.11)  through  (3.15)  we  see  that 


(3.14) 


(3.15) 


(3.16) 


With  this  T\  we  can  define  a  group  representation  on  X  as  the  direct  product  of  Tz  and 
T>;  that  is 

T(g)x  =  ( Tz(g)z,Tx(g)\)yx  e  X,Vg  €  G,  (3.17) 

which  then  implies  that 


F(T(g)x)  =  S(g)F(x),Vx  €  X,V<?  €  G.  (3.18) 


The  next  theorem  concerns  the  relation  between  the  original  representation  Tz  in 
(3.10)  and  the  induced  representation  T\  in  (3.16). 

Theorem  3.3  :  If  F(x)  =  0  and  z  =  Tz(g)z  for  some  g  &  G,  then  A  =  T\(g) A. 

The  computation  of  foldpoints  generally  involves  some  augmented  equation  which 
is  then  solved  by  a  Newton-like  method.  The  cutset-idea  introduced  in  [7]  provides  a 
natural  geometric  interpretation  of  augmentations.  In  particular,  in  [7]  it  was  shown  that 
with  certain  cuts  of  the  solution  manifold  it  is  possible  to  recover  many  of  the  results  of 
bifurcation  theory  required  in  computational  considerations. 

In  general,  we  consider  cuts  of  the  solution  manifold  M  by  the  orthogonal  complement 
IbJ-  of  any  given  subspace  W  of  X.  Thus  the  cut  orthogonal  to  W  through  the  given  point 
xo  €  M  is  the  intersection  M  D  (xo  +  W^)  and  therefore  consists  of  the  regular  points 
x  6  D  which  satisfy  the  augmented  equations 


F(x)  =  0,  IIw(x  -  x0)  =  0, 


(3.19) 


where  II w  is  the  orthogonal  projection  of  X  onto  W.  If  we  use  the  tangential  local- 
coordinate  system  (2.6)  of  M  at  x0  and  choose  W  as  iVo  =  An NXoM  of  (2.8),  then  (3.19) 
is  equivalent  with 

UNoDF(xoyr,(t)  =  0,t€  J,  (3.20) 

and  this  suggests  the  definition  of  the  cut  function 

/  :  J  C  TXoM  ~  No 
f(t)  =  UNoDF(xoYrj(t),  t  6  J. 

The  following  results  address  the  group  equivariance  of  the  cut  function. 


(3.21) 
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Lemma  3.4  :  For  any  point  x0  €  XG,  and  with  N0  =  A  D  NXoM,  we  have  T(g)(N0 )  C 
No,Vg<=G. 

The  proof  is  self-evident  and  need  not  be  detailed  here. 

Theorem  3.5  :  If  No  =  A  n  NXoM  {$}  and  IIjv0  denotes  the  orthogonal  projection 
onto  N0,  then  T(g) UNo  =  nNoT(g),Vg  €  G. 

Proof  :  Since  X  =  Z  ®  A,  and  X  =  NXoM  ®  TXoM,  it  follows  that 

I  =  2©(AnTIOM)©(AnlVIOM).  (3.22) 

Let  IV2  =  Z  ®  (A  fl  Tro  M),  Wi  =  A  fl  NXqM  =  N0.  Then  Lemma  3.4  implies  that  W}  is 
invariant  under  T  and  because  of  Z  =  A±,NXoM  =  (T^M)-1-  we  have  W2  =  W^.  Hence, 
W2  is  also  invariant  under  T  and  (2.22)  of  Theorem  2.5  applies  to  P  =  Un0.  | 

Theorem  3.6  :  Let  xo  €  M  fl  XG .  The  cut  function  f  defined  hy  (3.21)  at  xq  satisfies 

f(T(g)t)  =  T(g)f(t)Ng  g  G,Vt  €  J.  (3.23) 

Proof  :  (a)  We  show  first  that  the  function  77  in  the  characterization  (2.6)  of  a  local 
coordinate  system  satisfies 

v(T(g)t)  =  S(g)rj(t)yg  e  G,Vt  G  J.  (3.24) 

In  the  proof  of  Theorem  2.2  ([19,  Theorem  4.4])  the  inverse  function  theorem  ensures  the 
existence  of  an  open  ball  J  =  B( 0,  (3)  with  radius  [3  >  0,  of  a  neighborhood  V  C  Y  of  the 
origin,  and  of  a  unique  Cr-function  rj  :  J  »— ►  V  such  that 

H(t,  7?(t))  =  F(x0  +  t  +  DF(xoyr)(t))  =  0,  Vt  G  J.  (3.25) 

Since  T  is  unitary,  we  have  T(g)J  C  J ,  and 

S{g)H{i,  i i(t))  =  S(g)F{x0  +  t  +  DF{xoyT](t)) 

=  F(x0  +  T(g)t  +  DF(xoyS(g)T](t ))  (3.26) 

=  H(T(g)t,S(g)ri(t))  =  0. 

On  the  other  hand,  by  replacing  t  G  J  by  T(g)t  G  J  in  (3.25)  it  follows  that 

=  0.  (3.27) 

Since  T(g)t  G  J ,  the  solution  of  (3.25)  for  T(g)t  is  unique  and  hence 

r}(T(g)t )  =  S(g)r}(t). 

(b)  From  (a)  and  Theorem  3.5,  we  obtain  that 

T(g)f(t)  =  T(g)UNoDF(xQyV(t) 

=  JlN0DF(xoyS(g)g(t)  (3.2S) 

=  UNoDF(x0)tr1(T(g)t)  =  f(T(g)t), 
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which  proves  the  theorem.  | 

It  was  shown  in  [7]  that  the  cut  function  /  always  satisfies 

/(0)  =  0,  Df(0)  =  0.  (3.29) 

Accordingly,  the  second  singularity  index  of  the  foldpoint  x0  €  M  is  defined  (see  [7])  as 
the  number  k  >  1  such  that 

Djf( 0)  =  0,  j  =  0, 1, . . . ,  k,  Dk+1f( 0)  ?  0,  (3.30) 

if  there  exists  such  a  k,  else  k  =  oo.  With  this,  the  triple  ( p ,  q,  k )  consisting  of  the  Fredholm 
index  p  and  the  two  singularity  indices  q  and  k  forms  a  classification  for  the  foldpoints  of 
M. 

For  the  determination  of  the  form  of  the  cuts  the  generalized  Morse  Theorem  due  to 
Buchner,  Marsden  and  Schecter  [1]  is  applied  in  [7].  There  it  was  shown  that  the  form 
of  the  cut  orthogonal  to  No  of  (2.8)  at  the  foldpoint  Xq  of  type  (p,q,k)  is  essentially 
determined  by  the  zeros  of  the  (k  +  l)-form 

Q  ■  J  C  TX0M  JV0, 

Qm  =  — ^ e  j.  (3-31) 

More  specifically,  if  Q  is  regular  on  its  zero  set;  in  other  words,  if  DQ(t )  is  surjective  for 
each  nonzero  t  €  <5_1(0),  then  locally  near  xq  the  cutset  is  C1-diffeomorphic  to  the  zero 
set  of  Q.  In  line  with  this,  a  foldpoint  of  M  of  type  ( p,q,k )  will  be  called  nondegenerate 
if  its  associated  ( k  +  l)-form  (3.31)  is  regular  on  its  zero  set. 

Theorem  3.7  :  The  ( k  +  1  )-form  Q  associated  with  the  cut  function  f  in  Theorem  3.6 
is  equivariant  under  G, 


Q(T(g)t)  =  T(g)Q(t),Vg  €  G,Vt  €  J.  (3.32) 

Proof :  By  Taylor’s  theorem  we  can  write 

/(f)  =  Q(<)  +  <’(ll*ll‘+1).  (3.33) 

and  for  t  ^  0,  we  have 

=  Q(T(g)t)  +  o(||r(9)t||*+I) 

=  <?(7W)  +  °(l|f||*+'), 

T(g)f{t)  =  T(g)Q(t)  +  oflllH*-"  )■  (3.35) 

Now  divide  both  sides  of  (3.34)  and  (3.35)  by  and  set  v  =  Then  for  t  —*  0  it 

follows,  in  view  of  (3.23),  that 

Q(T(g)v )  =  T(g)Q(v),Vg  €  G,v  e  TXoM, \\v\\  =  1,  (3.36) 
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and  by  homogeneity  the  condition  ||u||  =  1  can  be  dropped  in  (3.36).| 

Remark  :  The  linearly  independent  vectors  v  €  Q~1( 0),u  ^  0  are  called  bifurcation 
directions.  By  Theorem  3.7,  if  v  is  the  zero  of  Q,  so  is  the  T(g)v,g  6  G,  that  means,  the 
representation  T  acts  on  Q-1(0)  by  permuting  the  bifurcation  directions  under  symmetry. 
Therefore,  linearly  independent  bifurcation  directions  v  £  Q-1(0),  ||u||  =  1  correspond  to 
different  orbits  through  x0  on  M . 

If  x0  is  afoldpoint  in  Xc DM  it  may  occur  that  the  bifurcation  directions  at  x0  belong 
to  XH  for  some  isotropy  subgroup  H  of  G.  In  that  case  the  symmetry  group  of  equation 
(2.1)  remains  unchanged  but  the  bifurcation  directions  spontaneously  break  symmetry  (in 
the  absence  of  an  external  symmetry  breaking  perturbation). 

An  important  question  in  the  study  of  symmetry  breaking  is  the  determination  of  the 
subgroup  H  of  G  which  breaks  the  G-symmetry.  For  this,  we  present  a  manifold  version  of 
the  equivariant  branching  lemma  which  provides  some  choices  for  the  symmetry  breaking 
subgroup. 

Theorem  3.8  :  Suppose  that  Xo  G  M  fl  XG  is  a  foldpoint  of  type  (p,q,  1),  and  H  is  an 
isotropy  subgroup  of  G.  Let  Z0  =  ZC\TXoM.  Ifdim(Z0nXH)  =  1,  say,  Z0 C\XH  =  span(u0), 
then  vq  is  a  bifurcation  direction  in  M . 

Proof  :  It  was  shown  in  [4]  that  if  x0  is  a  foldpoint  of  type  ( p ,  1, 1);  that  is,  if  dimZo  =  1, 
and  Z0  =  span(uo),  then  u0  is  the  bifurcation  direction.  The  key  point  here  is  to  reduce 
the  original  multiple  foldpoint  problem  of  type  (p,q,  1)  to  the  Tf-reduced  simple  foldpoint 
problem  of  type  (pn,  1,1).  If  we  can  show  that  Qh(v0)  =  0  then  by  the  equivariance  of  Q 
(Theorem  3.7),  we  also  have  Q(v o)  =  0.  More  specifically,  we  need  to  show  that 

qH  =  dim  ( TZ0M„  n  ZH)  =  1,  (3.37) 

where 

ZH  =  {z  €  Z;  T(h)z  =  z,Vh  6  H}  (3.38) 

is  the  fixed-point  subspace  of  H  on  Z.  Since,  by  definition,  Fh{x )  =  F(Phx),  it  follows 
that 

DF„(x)  =  DF(PHx)PH  =  DF(x)PHyx  €  XH ,  (3.39) 

which  proves  the  identity 

ZH  n  ker  DF„(x)  =  XH  D  (Z  fl  ker  DF{x))  =  XH  D  Z„.  (3.40) 

This,  together  with  the  assumptions  xo  G  M  fl  ArG  C  M  fl  XH  and  dim(Xw  fl  Zo)  =  1, 
implies  that  (3.37)  holds.  | 

Remark  :  It  can  be  shown  that  if  Q  is  regular  at  Vo,  then  the  bifurcated  branch  tangent 
to  vo  preserves  the  symmetry  of  the  isotropy  subgroup  of  i?o  (see  [16]). 

For  the  numerical  implementation  of  Theorem  3.8  note  that 

Zo  =  Zf]TXoM  =  ker  DzF(x0),  (3.41) 
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and 


XH  nZ0  =  z”  =  {z  €  Z0;  T(h)|2o*  =  z,Vh  6  #}.  (3.42) 


If  (Pq)h  denotes  the  orthogonal  projection  from  Z0  onto  Zq  then  by  (2.18)  and  (2.19)  we 
have 


dim(Z0")  =  _L  Y,  MIWU) 

^  '  h£H 

=  tr  (7WT  S  =  tr(Ci’o)//)- 

1  1  heH 


(3.43) 


Lemma  3.9  :  (see  Chapter  I, §2  of  [15])  Let  A  be  a  linear  operator  on  a  finite-dimensional 
space  U  and  V  a  subspace  ofU  that  is  invariant  under  A.  Moreover,  let  P  be  a  projection 
of  U  onto  V  and  Ay  the  restriction  of  A  to  V.  Then  we  have 


tr(Ay)  =  tr  (AP). 


(3.44) 


Let  ( Pz)h  be  the  orthogonal  projection  from  Z  onto  Zh  and  11^,  DtF(x0)  the  orthog¬ 
onal  projection  from  Z  onto  ker  DzF(xq).  Then,  by  Lemma  3.9,  (3.43)  becomes 


dim(A'H  n  Z0)  =  tr((Pz)HflketDzF(x0)),  (3.45) 

where 

(Pz)h  =  E  <3-46) 

h£H 

and  IIker  dzF(x0)  *s  available  from  the  singular  value  decomposition  (SVD)  of  DzF(x o)  in 
the  foldpoint  computation  (see  later). 

To  determine  which  isotropy  subgroup  H  satisfies  the  condition  dim(Z0  PI  XH)  =  1, 
we  usually  search  through  the  lattice  of  the  isotropy  subgroups  of  G  starting  from  the 
root  G.  Note  that  when  dimA-77  =  1  and  XH  D  Zq  ^  {<j>}  then  qu  =  1.  An  isotropy 
subgroup  H  is  called  maximal  if  it  is  not  a  proper  subgroup  of  any  other  isotropy  subgroup 
except  G.  If  dim.Yw  =  1  then  H  is  a  maximal  isotropy  subgroup.  A  major  advantage 
of  considering  such  maximal  isotropy  subgroups  is  that  we  can  use  (2.18)  to  determine  all 
maximal  isotropy  subgroups  H  with  dim  XH  =  1  and  then  use  them  as  our  first  choice  to 
test  whether  Theorem  3.8  is  satisfied  for  H  (see  Sattinger  [22]  and  Golubitsky  ei  al  [11]). 

In  view  of  Theorem  3.8,  Golubitsky  [9]  conjectured  that  generically  the  solutions  of 
symmetry  breaking  bifurcation  problems  all  satisfy  dim(ArW  D  Zq)  =  1.  It  is  now  known 
that  the  conjecture  is  false.  Counter-examples  for  higher  dimensional  representations  of 
50(3)  have  been  found  by  Chossat  [2]  and  for  higher  dimensional  representations  of  0(3) 
by  Lauterbach  [13].  On  the  other  hand,  Field  and  Richardson  [8]  show  that  for  finite 
groups  generated  by  reflections  the  conjecture  is  true. 

The  above  discussion  about  finding  the  symmetry  breaking  subgroup  H  of  G  is 
strongly  representation  dependent.  If  we  change  the  representation  we  will  need  to  rebuild 
the  list  of  choices  for  H .  Thus  the  question  arises  whether  there  exists  any  common  list  for 
different  representations.  A  positive  answer  has  been  given  by  Rabier  in  [17]  and  [18].  More 
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specifically,  it  was  shown  that  Lie  groups  possess  subgroups  of  a  special  type,  called  intrin¬ 
sic  isotropy  subgroups,  which  are  characterized  by  a  representation-independent  property 
provided  the  representations  are  assumed  to  admit  covariant  linear  isomorphisms  with  neg¬ 
ative  determinant.  Such  intrinsic  isotropy  subgroups  play  an  important  role  in  symmetry 
breaking  bifurcations. 

Let  T  be  any  Lie  group  and  T'  a  subgroup  of  T.  Then  T'  is  an  intrinsic  isotropy 
subgroup  if  for  any  representation  T  of  T  in  C?jL(Rn)  with  given  n  >  1,  any  T-covariant 
linear  isomorphism  L  G  GL(Rn)  satisfies 

sgn  det  Lp/  =  sgn  detZ.  (3.47) 

Suppose  that  the  Lie  group  T  is  given  along  with  one  of  its  representations  T  in 
GZ(R").  We  consider  a  C1  mapping  g(=  g(p,x))  :  R  x  Rn  h- ►  Rn  such  that  g(p,  0)  =  0 
and  assume  that,  for  every  /x  £  R ,g(p,*)  is  T-covariant.  Moreover,  we  assume  that 
det  Dxg(p,  0)  changes  sign  as  p  crosses  0.  With  these  hypotheses,  the  following  results 
hold. 

Theorem  3.10  :  ([17,  Theorem  2.1])  Let  T'  C  T  be  a  (maximal)  intrinsic  isotropy  sub¬ 
group.  Then,  there  is  a  branch  of  nontrivial  solutions  of  g(/x,  x)  =  0  with  T'-symmetry 
that  bifurcates  from  (0,0). 

For  the  existence  of  nontrivial  intrinsic  isotropy  subgroups  (the  trivial  case  is  T'  = 
{/}),  it  was  shown  in  [17]  that  a  compact  Lie  group  T  has  a  nontrivial  intrinsic  isotropy 
subgroup  if  and  only  if  it  is  not  isomorphic  to  a  direct  product  of  copies  of  Z 2.  If  the 
elements  of  odd  order  in  T  form  a  subgroup  T0ad  (e.g.  if  T  is  supersolvable),  then  the 
subgroup  (T2)  is  an  intrinsic  isotropy  subgroup  of  T,  where  (T2)  is  the  group  generated  by 
elements  of  the  form  72,7  €  T.  Surprisingly  the  converse  is  also  true  (see  [IS]). 

Now  we  are  ready  to  adapt  the  computational  method  for  (p,  1,1)  foldpoints  to  the 
reduced  case  of  (p/f,  1, 1)  foldpoints.  For  this,  we  need  the  additional  assumption 

Assumption  3.11  :  There  exist  a  point  xq  €  M  H  XG  which  is  not  a  foldpoint  with 
respect  to  A. 

Lemma  3.12  :  Under  the  Assumption  3.11  we  have 

dim  ZH  =  dimFw,  and  dimker  DzFh(x)  =  dimker  DzFnix)1.  (3.48) 

Let  x*  G  XHflM  be  a  non-degenerate  foldpoint  of  M  with  respect  to  A  of  type  (p,  q,  1) 
where  q  >  1.  Then  we  have  rank  DzF(x*)  =  m  —  q  and  hence  rank  DzF(x )  >  m  —  q 
for  all  points  x  in  some  open  subset  Uq  of  X  containing  x*.  Suppose  further  that  qu  — 
dimker  DzFh(x*)  =  1  and  set  Uq1  =  Uo  D  XH .  Then  we  have  rank  DzFh(x*)  =  mh  —  1 
and  hence  rank  DzFh(x)  >  mu  —  1  for  all  points  x  in  some  open  subset  Vo  of  U q1  containing 
x*.  In  other  words,  for  any  foldpoint  x  G  Vo  DM  the  first  singularity  index  of  F  is  at  most 
q  but  the  first  singularity  index  of  Fh  equals  1. 

Let  Lh  be  an  orthonormal  matrix  whose  columns  span  ZH ,  then  from  (3.48)  the 
matrices  DFh(x*)Lh  and  ( DFh(x*)LhY  have  one-dimensional  nullspaces.  Let  u*H,c*H  G 
Rm  be  any  given  vectors  such  that 

CH  &  rge  DFh(x*)Lh,  u*h  &  rge  (DFh(x*)LhY ,  (3.49) 
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and  for  any  x  €  Vo  consider  the  matrix 


Ah{x) 


(  DFh(x)Lh  c*h\ 

V  («*)*  o  ) ' 


(3.50) 


Then  there  exists  an  open  subset  Vj  of  Vo  containing  x*  such  that  Ah{x )  is  non-singular 
for  all  x  €  V\.  Now  the  solutions  of  the  two  linear  systems 


Ah(x) 


u°H(x) 

Vh(x) 


=  0,  Ah( 


\7H(X) 


+ 


—  0,  i  €  Vj 


(3.51) 


are  uniquely  determined  for  each  x  G  V\.  Applying  Theorem  1  and  Theorem  3  of  [4]  to 
the  reduced  (ph,  1,1)  foldpoint  we  have 


Theorem  3.13  : 

of 


There  exists  an  open  neighborhood  C  V  of  x*  such  that  the  solutions 


(3.52) 


are  exactly  the  foldpoints  of  M  in  Mh  fl  Vj  with  type  (ph,  1,1)  and  these  solutions  form 
a  (pn  —  1) -dimensional  submanifold  of  Mu- 


Remark  :  Note  that  usually  pn  <  p,  if  pu  =  2,  then  the  foldpoint  submanifold  with 
symmetry  H  is  one-dimensional  and  we  can  apply  PITCON  (see  [19])  to  trace  this  sub¬ 
manifold.  On  the  other  hand,  if  pH  >  3,  we  can  apply  MATRIG  (see  [20])  to  compute  a 
simplicial  approximation  of  an  open  subset  of  this  sunmanfold. 


We  summarize  the  results  in  the  form  of  the  following  algorithm. 

Algorithm  3.14  :  Let  x  be  an  approximate  foldpoint  of  type  (p,q,  1)  of  M  with  respect 
to  A: 

(1)  Compute  the  singular  value  decomposition  (SVD)  of  DzF(x ) 


B(x)lDzF(x)A{x)  =  diag(cri(x), . . .  ,crm(x)). 


(2)  Let  a'(x),i  =  1  be  the  columns  of  the  matrix  A(x)  and  compute 

the  projection 


nicer  DtF(x)  —  AqAlq,  Aq  =  (am  ,+  1,...,am). 

(2)  For  all  possible  symmetry  breaking  isotropy  subgroups  H 
3a)  compute  the  projection 


(Pz]H  =  jh  £  T,W; 


hEH 


3b) 


and  determine  which  isotropy  subgroup  H  satisfies 
dim(XH  fl  ker  DzF(x))  =  tr((Pz)//nkerDzF(lo))  =  1. 


14 


(4)  Compute  axH  =  (Pz)ho.'(x)  for  2 
some  io,  then  take  u*H  =  a$. 

(5)  Compute  the  projection 


=  m, . . . ,  m  —  q  +  1  until  a ^  0  for 


Qx  = 


(6)  Let  bi(x),j  =  1 ,m  be  the  columns  of  the  matrix  B(x )  and  compute 
b]H  =  Qh^{x)  for  j  =  m, ...  ,m  —  q  +  1  until  bjf  ^  0  for  some  jo,  then  take 


C*H  ~  b}H- 


(7)  Set  xo  =  x  and  for  k  =  0,1, ... ,  until  convergence 
7a)  solve  the  systems 

«■»(;;£!)*(!) -*  ■*"<■•>•  (SSI)* 

7b)  compute  G h(xic)  and  DGn(xk) 

7c)  and  compute  the  next  iterate 

Xk+i  —  Xk  -  DGH{xk)-lGH{xk) 


4.  Numerical  example.  A  programming  package  BISYM  has  been  written  which  im¬ 
plements  the  methods  of  Section  3.  In  order  to  indicate  the  performance  of  the  package 
BISYM,  we  consider  the  four-box-Brusselator  reaction  problem  (see  [5])  in  unfolded  form 


where 


AAit  —  G(u,  v)  +  8id  +  62b  =  0 

lOXAv  -  H(u,  v)  +  O.S1533(^a  +  S2b)  =  0 


2-10-1 

-12-10 

0-12-1 

-10-12 


u,  v  6  ft4 


G(u,v)  =  {g(uuvi),...,g(u4,v4)y,  g(s,t)  =  2-6.9  s  +  s2* 
H(u,v)  =  {h{u\,v\ h(u4 ,  u.*))*,  h(s,t)  =  5.9  —  s2t 
a  =  (0,1, 0,-1)*,  6  =  (1,0,-1,0)*. 
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With  x  =  (u,v,  \, 61,82)  €  R"  and  n  =  11, m  =  8,p  =  3  this  is  a  problem  of  the  form 
considered  in  sections  2  and  3.  It  is  evident  that  the  equivariance  relation  (2.12)  is  satisfied 
with  a  group  G  that  is  isomorphic  to  the  dihedral  group  D4  : 


G  —  {I,  Ri,  R2,  R3,  Si,  S2,  S3,  S4} 

I,  Ri ,  R2 ,  R3 ,  Si ,  S2 ,  S3 ,  S4  are  permutation  matrices  representing  rotations  and  refections 
(see  Figure  1).  The  induced  group  representation  T\  with  respect  to  this  unfolding  is  given 
by 


Tx(Ri)  = 


A  0  0\  A  0  0  \  A  0  0\ 

0  0  1  TA(U2)=  0  -1  0  J  T\(R3)  =  (  0  -1  0 

\o  -1  0/  \o  0  -1/  \o  0  1/ 


and 


A  0 

T\(Si)  =  0  0 

\0  1 

A  0 

Ta(53)=  0  1 

\0  0 


0\  (1  0  0  \ 
1  Tx(S2)=  0  0  -1 

0/  \0  -1  0  / 

0  \  /I  0  0\ 

0  Tx(S4)=  0-10 

-1)  \ 0  0  1/ 


(4.4) 


(4.5) 


The  subgroups  of  G  are  : 


G  —  {I,Ri,  R2,  R3,  Si,S2,  S3,  S^^r  —  {I,Ri,R2,R3}i 

'Em  =  {I,R‘2,Si,S2},Ed  =  {I,  R2,  S3,  S4},T,P  =  {I,  R2},  (4-5) 

Si  =  {I,Si},E2  =  {/,52},E3  =  {I,S3},E4  =  {/,54},E  0  =  {/}. 


It  can  be  shown  that  only  G,  E4 ,  E2,  S3,  E4  and  Eo  are  isotropy  subgroups  and  :r0  = 
(u* ,  v* ,  A*,  0, 0)  with 


e  =  (1, 1, 1, 1)',  u*  =  2e,  v*  =  2.95e,  A*  =  0.04536 


(4.7) 


is  a  foldpoint  of  (4.1)  with  first  and  second  singularity  indices  2  and  1,  respectively.  By 
Theorem  3.8  and  the  formula  (3.45),  we  can  check  that  symmetries  breaks  from  G  to 
E i,i  =  1 , . . . ,  4  respectively  at  xq . 
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w  w 


Using  (4.4),  (4.5)  and  (2.18)  we  obtain  dim Mjv  =  2 ,i  =  1,...,4  (note  dim A/V,  = 
dim  AE' ).  Applying  the  subroutine  PITCON  of  BISYM  we  can  trace  the  foldlines  in 
symmetry  E,-,  *  =  1, . . . ,  4. 


Figure  2 


Figure  2  gives  the  projections  of  the  foldlines;  that  is,  the  caustics  of  the  symmetries 
Ei,  E2,  S3,  and  E4  in  the  parameter  space  (62, 61,  A).  Each  caustic  has  the  shape  of  a  cusp. 
It  may  be  interesting  to  note  that  the  opening  of  the  cusp  of  Si,  E2  is  opposite  to  that  of 


(a)  (b)  (c) 


Figure  3 

Figure  3  provides  pictures  of  the  projections  of  the  foldlines  of  the  symmetries  E2  and 
4  onto  various  2-dimensional  spaces,  they  indicate  that  the  projections  of  foldlines  of 
2  and  E4  onto  the  space  of  (vi,  A,  62)  behave  as  the  cusp  functions  v\  +  A  *  V\  4-  ^2  and 
t’j  —  A  *  v\  + 62 ,  respectively. 
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Applying  the  subroutine  MATRIG  of  BISYM  we  can  also  compute  the  simplicial  trian¬ 
gulation  of  the  reduced  manifold  A/v^ , . . . ,  Ms4 ,  respectively,  in  the  coordinates  (i>j ,  A,  ) 
(see  Figure  4). 


Figure  4 


Figure  4  certainly  indicates  how  complicated  the  solution  behavior  is  in  this  example. 

Acknowledgment.  I  would  like  to  express  my  appreciation  to  my  advisor,  Professor 
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